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This report attempts to capture the essential workings of gauge links 
(Wilson lines) inside gauge-invariant formulations of parton distribution 
functions in QCD and gain some deeper insight into their key (renormal- 
ization) properties. We show, in particular, that the one-loop anomalous 
dimension of the Cherednikov-Stefanis quark TMD PDF is in the lightcone 
gauge A + = 0, combined with the Mandelstam-Leibbrandt pole prescrip- 
tion, the same as that obtained in the special covariant gauge a = —3, 
leaving no uncanceled rapidity singularities. 

PACS numbers: ll.15.-q, ll.10.Gh, 12.38.Aw, 12.39. St 

1. Introduction 

The standard way to remove the gauge-dependence of nonlocal corre- 
lators in gauge theories, like QCD, is to include path-ordered exponential 
factors of the gauge field that are absent in the original Lagrangian of a local 
quantum field theory. These operators are in general contour-dependent and 
give rise to new divergences, called in modern jargon rapidity divergences, 
that are not related to the standard singularities — ultraviolet (UV) and 
infrared (IR) — that appear in Feynman diagrams. They originate from con- 
tour irregularities caused by topological obstructions — endpoints, cusps, 
self-crossings — of the gauge contours entering the exponents of the gauge 
links and affect the renormalization properties of the QCD correlators pQ. 
The trouble is, such divergences have to be regularized and there is great 
ambiguity in adopting an appropriate subtraction procedure within a valid 
factorization scheme. Moreover, the use of the lightcone gauge quantization 
depends crucially on the adopted boundary conditions imposed on the gluon 
propagator in order to treat the gauge links at infinity. 
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These issues will be considered here in some detail. In Sec. HJ we will 
address the one-loop virtual corrections of the quark propagator in a general 
covariant gauge within two different gauge- invariant schemes: The Mandel- 
stam formalism [2] and the z-field algorithm O 0]. Then, in Sec. [3j we 
will discuss — as an example of a transverse-momentum dependent (TMD) 
parton distribution function (PDF) — the quark in a quark TMD PDF, 
fq/q-, using the lightcone gauge subject to various boundary conditions on 
the gluon propagator. Finally, the conclusions will be presented in Sec. HI 



In this section, we consider a gauge-invariant mesonic type correlator in 
QCD, viz., 



is a path-ordered gauge link along some arbitrary contour V between x\ 
and X2 and tfj = (1/2)A^ are the Gell-Mann matrices of SU(3) C - In prin- 
ciple, any contour between the points x\ and X2 is admissible. Contour 
obstructions give rise to rapidity divergences and hence contribute to the 
anomalous dimension of the correlator. For our considerations in this sec- 
tion, we assume that the contour is smooth. Because only the two endpoints 
entail rapidity divergences, it is sufficient to employ the straight line T join- 
ing x\ and x%. The reason is that other features of smooth contours, e.g., 
their length, do not affect the renormalization properties of the twist-two 
mesonic correlator and hence f can be used as a subtraction contour for its 
renormalization for all members of the universality class of smooth contours. 
Differences in the definition of the mesonic correlator for different smooth 
contours show up at the next higher twist level [5]0 

The gauge link can be treated in two different ways. One can evalu- 
ate the path-ordered exponential as a power series in the coupling using 
the Mandelstam formalism [2J. This approach [B] will be on focus in the 
first subsection below. Another option is to apply the z-field algorithm 
[31 H], which is based on an effective Lagrangian describing the interaction 
of one-dimensional auxiliary Fermion fields with the gluon field and trade 
contours for "particle trajectories". Fully quantized results are finally ob- 
tained by performing a functional integral over the z-field fluctuations to get 
the mesonic correlator in second quantization. Our presentation here follows 
the analysis of [7] with more details to be given in a future publication. 

1 I thank Sergey Mikhailov for discussions on this point. 



2. Mesonic correlator 
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Fig. 1. One-loop contributions to the vacuum expectation value of the mesonic 
correlator in Eq. ([3]) within the Mandelstam formalism, (a) Usual quark self-energy 
with gluon attachment denoted by a curly line; (b) and (c) contributions due to Eqs. 
((3]) and ([3]), respectively. Curved solid lines denote quark propagators, whereas the 
gauge contours are represented by dashed straight lines with the symbol ® denoting 
line integrals along them. 



2.1. Mandelstam formalism 

We carry out the vacuum expectation value of M(x%, X2\T) to the second 
order of the unrenormalized coupling constant g, i.e., 

(0|T (q(x 2 )[x 2 , xtlTMx^lO)® = (0\T (q(x 2 )[x 2 , x^qix^O}^ 

+ (0|T(g(x 2 )[x 2 ,x 1 |r]( 1 )g(x 1 ))|0)« 

+ (0|T(g(x 2 )[x 2 ,x 1 |r]( 2 )g(x 1 ))|0)(°)(3) 

with each contribution evaluated to the appropriate order. The first term 
is the usual gauge-dependent quark self-energy to 0(g 2 ), whereas the sec- 
ond and the third term represent the 0(g 1 ) and the 0(g 2 ) contributions 
stemming from the gauge link (termed in [6] the connector), respectively: 

[x 2 ,xi|r]W = ig / dZpA^z) , (4) 



(inV f X2 f X2 

lx 2 , Xl \T}W = m_ / dz l d z' v V{A^z)A v {z!)) . (5) 
21 Jxi[r] Jxi[t] 

The three terms displayed in Eq. Q are shown graphically in Fig. [TJ 

It was shown in |6] that the rapidity divergences entailed by the end- 
points in diagrams (b) and (c) in Fig. [1] can be controlled by dimensional 
regularization with no need to involve additional regulators. The explicit 
calculation can be found there. Here we only quote the final results recalling 
that we are dealing with quark fields that are Heisenberg field operators in 
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the interaction picture (where they are free operators) so that one has to 
perform all Wick contractions in 



M(x 1 ,x 2 \T) = ( 



T^q (x 2 )[x2,x 1 \T}qo(xi)exp (i J d 4 y£ int (y) 



.(6) 



The mesonic correlator in a general covariant gauge a can be written in the 
form 

M^ j \p\T) = M^' b \ P ) + (1 - a)M ( 2 a ' b) {p) (7) 

with 

D F(k) = -J-^ + (l-a)^-^, S F (P) = . (8) 
k z + le k z + ze p — m + te 

where the first superscript in Ad^'i' indicates the order of the expansion of 
the gauge link, whereas the second one denotes the order of the coupling 
taken into account in the Wick contractions. One finds [6J that the connector 

parts are interrelated: M%~ (p) = —2M^^(p), while the remaining part 

of M$ (p) cancels the (1 - a) contribution of M.f ,2 \p), so that M^ 2 \p\T) 
is gauge- parameter independent and multiplicatively renormalizable: 

M(p\f)U = ^h y teid^(p|f)lbare ■ (9) 

Using dimensional regularization we obtain in the MS scheme the fol- 
lowing renormalization constants (D = 4 — e) [B] 

^hybrid = Z 2q Z con = 1 + |^<? F i + 0(g 4 ) (10) 

Z con = 1 + #^ F (3 + a)- + 0{g A ) , (11) 

Z 2q = l-|^C F ai + 0( 9 4 ) , (12) 

where Cp = (JV? — 1)/(2JV C ) = 4/3 for N c = 3 The associated anomalous 
l,,d_ 



dimensions (7 = \p-§-Yn.Z) are given by 



7hybrid = l2q + 7con = ~J^2 C F + Otf) ( 13 ) 



,2 



7con = ~J^2 C ^ + «) + , (14) 



-,2 



729 = T^2 CFa + °( 9 > (15) 

and bear no dependence on the geometric features of the contour, e.g., its 
length or its derivatives d n x{r) / '(dr) n . 
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Fig. 2. One-loop contributions to the mesonic correlator in the .z-field formalism: 
(a) z-field (broken line) self-energy; (b) and (c) vertex corrections from zA^^^p-z. 



2.2. z-field formalism 

Though the correlator M{x\,X2\F) is multiplicatively renormalizable, it 
cannot be written contour-independently in the factorized form M(x±, x 2 \T) = 
Q(x2\T2)Q(xi\Ti) starting from the QCD Lagrangian in terms of the Man- 

delstam field Q(xi\Ti) = Pexp ig /*J r ri i dz^Aa(z)t a q{x\). However, for 

smooth contours one can factorize the connector [xi,X2|T] according to the 
algebraic identity [B,A\T] = [B,C\T 2 ][C,A\T 1 ] with r x U T 2 = T and then 
shift C — > oo. This trick allows factorization and multiplicative renormal- 

— 1/2 

ization, i.e., Q(x|r)| ren = Z hybrid Q(x|Ij|bare- A more far-sighted option is 
to trade gauge contours in favor of "trajectories" of fictitious particles de- 
scribed by the following effective Lagrangian 



eff 



da 



z(a)-^z(a) 



igz(a)^^A»(x(a))z(a) 



iXz + izX 



(16) 



and supplement £qcd by two additional terms pertaining to two extra Feyn- 
man rules: one for the z-field propagator and the other for the z-field-gluon 
vertex (see [7]). This allows one to write the gauge link as a path integral: 



V exp 



1 9 



CT2 J r 

daA£(x(a))t a ^ 

01 



da 



(0|^(a 2 )%(ai)|0) 



(17) 



f° 2 daz(a)D a z(a) . What is more, perform- 



where (Q) = J VzDzQexp (rr 

ing now the calculation of the radiative corrections to the mesonic corre- 
lator within this approach (see Fig. [2]), one finds at 0(g 2 ) that the local 
combination z(a x )q(x) (analogous to the nonlocal field Q(x\T) in the Man- 
delstam formalism) gets renormalized by the renormalization constant [7] 



Z 



-1/2 



z 



-1/2 7 -l/2 



Z*t", where Z 3z 



1 + |^C F (3 



a)- resembles Z c 



1/2 

CF ~~ ^CF ^2q ^3z ' 

yielding to the same anomalous dimensions as before. As a result, the 
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Slavnov- Taylor identities to 0(g 2 ) are fulfilled (Ca = N c = 3): 

Z con Zi q Z 3 Z 3 Z- iz C^g 2 , . 1 . . 

- 1 + ^r( 3 + «)-- ( 18 ) 



J hybrid 



Z x Zi Z l2 32vr 2 



e 



One can now use the z-field formalism and perform a short-distance expan- 
sion of M(x\, x 2 \T) for x\ — > x 2 : 

M(x 1 ,x 2 \T) = q(x 2 )[x 2 ,x 1 \T]q(x 1 ) = q(x 2 )z(a 2 )z(a 1 )q(xi) (19) 

« ^cW^^^^.-.^oS.^^c) , 

where (cc2 — xi) 2 = z 2 ~ 0, 7^ and the composite non-singlet quark 
operator of lowest twist reads 

^ y ^ y 

' zn ...z^n v*"i w ; ~~ "A-v /mi mi ■ ■ ■ mjv 1 



o2 1 ..^ w (a!,C) = g(a!)7 w ^ 1 ...^[x > x|C] 9 ( a :) . (20) 



Here [x, x|C] is a gauge link along the closed loop C = TUT' and the short- 
distance expansion is valid for any point x because the smooth contour T' 
can be stretched to 00 by virtue of the independence of the renormaliza- 
tion constants on L(T'). An immediate important conclusion is that in the 
special gauge a = —3 all contour- (or z-field-) related divergences cancel 
among themselves so that the residual renormalization effects can be ab- 
sorbed into Z3, while the Ward identity Z\ q = Z 2q is preserved like in QED. 
It is expected that going to the next higher loop, one will obtain a similar 
result but for a* = —3 + 0(g 2 ) as found in [8] in the context of multiloop 
contributions to the nonsinglet QCD evolution equations. 



3. Gauge-invariant correlators for TMD PDFs 

With such issues in mind, let us focus attention on the question of the 
appropriate definition of a TMD PDF (e.g., [9] and references cited therein). 
In [10] we have shown that f q / q (x,k±) given in [11] cannot be regularized 
completely using dimensional regularization in the lightcone gauge A + = 
(A-n~) = with (n~) 2 = in conjunction with the retarded (ret), advanced 
(adv), or principal- value (PV) pole prescription on the gluon propagator: 

^VW" q 2_ X 2 + iQ \9?» [g+] I > W 



where 
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The reason is that the residue of the e pole contains a rapidity divergence 
that entails an extra anomalous dimension and calls for an additional sub- 
traction procedure. This can be achieved in terms of a soft factor R that 
has to be included into the definition of the TMD PDF [ID]. Its anomalous 
dimension 7^ serves to cancel £7 with the effect that the anomalous dimen- 
sion of fg, coincides with the result one would obtain in a covariant gauge 
for a direct smooth contour between the two field points i.e., Eq. (13). It 
turns out [ID] that this anomalous-dimension artifact has at one loop the 
same structure as the universal anomalous dimension of a cusped contour 
[12] . 5~f = — ^Cphi^r, which becomes infinite when 77 — >■ 0. For such 

^^-independent pole prescriptions the gluon propagator is not transverse: 
n^D^ 7^ 0. On the other hand, it was shown in [13] that using instead the 
(/^-dependent Mandelstam-Leibbrandt (ML) [14] pole prescription 



one has n^D^ L = 0. The opposite claims by Collins in [15] are likely the 
result of misunderstanding, as his equation (15) shows the gluon propaga- 
tor subject to the Principal Value pole prescription, i.e., Eq. (|22j) . This 
gluon propagator is indeed not transverse. The same applies to those 
obtained with the advanced or retarded pole prescription, because in all 
these cases the transverse gauge field A J -(oo~;^j_) = £^C oc 'V- L In A|£j_| is 
not purely transverse (in contrast to the ML case) but depends through 
Coo = {0(adv), — l(ret), — i(PV)} on the imposed boundary conditions. 

Another issue raised by Collins in [15] is whether the definition of f^j 
contains uncanceled divergences originating from the self-energy of the gauge 
links. We have shown in |10t [T3] that all UV divergences from the q mo- 
mentum integrations can be regularized dimensionally and give 1/e poles, 
while collinear poles are controlled by the quark virtuality p 2 < 0, and IR 
singularities are regularized by an auxiliary gluon mass A that drops out at 
the end. On the other hand, overlapping divergences have, in general, to be 
cured by a subtraction procedure encoded in the soft renormalization factor 
R and appear in In ^ in terms of an auxiliary mass rj. In the A + = gauge 
all diagrams with gluon attachments to the gauge links (longitudinal or 
transverse) either vanish identically or cancel partly against contributions 
from cross-talk diagrams with gluon attachments between the quark line 
and a gauge link. For the adv, ret, and PV pole prescriptions, all remain- 
ing divergences are taken care of by the soft factor rendering f^, regular. 
Employing the ML prescription, one gets a result that is reminiscent of 
that obtained in the special covariant gauge a = —3 in which all rapidity 




(23) 
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divergences emerging from the endpoints of smooth contours cancel among 
themselves. At higher loops, the gauge a = —3 will receive a s corrections 
with coefficients that can be determined by demanding the validity of the 
QED-like Ward identity Z\ q = Zi q . 

4. Conclusions 

We have shown the multiplicative renormalization of the gauge-invariant 
mesonic vacuum correlator in the nonlocal Mandelstam approach for smooth 
contours and found that to 0(g 2 ) it is equivalent to the result obtained 
in the local z-field effective formalism. At one loop, all contour- (or z- 
field) related rapidity divergences cancel among themselves using the special 
gauge a = —3. We argued that a proper definition of the quark TMD 
PDF must account for the appropriate subtraction of a rapidity divergence 
that overlaps with the usual UV singularities and cannot be regularized 
dimensionally. Employing the lightcone gauge with (/"-independent pole 
prescriptions, this can be achieved via a soft renormalization factor [10j . 
The imposition of the (/"-dependent Mandelstam-Leibbrandt prescription 
removes all rapidity divergences and reproduces at one loop the results of 
the special covariant gauge a = — 3 with the soft factor reducing to unity. 
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